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Discovery of the Youngest Molecular Outflow associated with an 
Intermediate-mass protostellar Core, MMS-6/OMC-3. 

Satoko Takahashi, 1 and Paul T.P. Ho 1,2 
ABSTRACT 

We present sub-arcsecond resolution HCN (4-3) and CO (3-2) observa- 
tions made with the Submillimeter Array (SMAQ|), toward an extremely young 
intermediate-mass protostellar core, MMS 6-main, located in the Orion Molecu- 
lar Cloud 3 region (OMC-3). We have successfully imaged a compact molecular 
outflow lobe (^1500 AU) associated with MMS6-main, which is also the smallest 
molecular outflow ever found in the intermediate-mass protostellar cores. The 
dynamical time scale of this outflow is estimated to be <100 yr. The line width 
dramatically increases downstream at the end of the molecular outflow (Av ~25 
km s _1 ), and clearly shows the bow-shock type velocity structure. The estimated 
outflow mass (f^lO" 4 M ) and outflow size are approximately 2-4 orders and 1-3 
orders of magnitude smaller, while the outflow force («10 -4 M km s _1 yr -1 ) is 
similar, as compared to the other molecular outflows studied in OMC-2/3. These 
results show that MMS 6-main is a protostellar core at the earliest evolutionary 
stage, most likely shortly after the 2nd core formation. 
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INTRODUCTION 



Molecular outflows are ub i quitously observ ed in low- to high- mass star forming regions 



(e.g.. 



Bachiller fc Taf allaj 1 1 9 9 9l ; lArce et al.ll2007l ). They have a variety of characteristics such 
as size (0.0 1 pc-a few pc) , mass (10" 4 -10 3 M ), and outflow force (KT 7 -! M km s" 1 



yr 1 ) (e.g., IWu et al.l I2004J ). These outflows are known to result from the entrainment of 



circumstellar gas, swept-up by the primary jet. 

It is well accepted that primary jets are launched by m agnetohydrodynamic (MHD) pro 



cesse s in the accretion disks surrounding a protostar (e.g., iPudritz et al.l 120071 : IShang et al. 



20071 ). They play an important role to remove angular momentum from the circumstellar 
disk and thereb y promote mass accre tion onto the central star. Furthermore, recent 3D-MHD 
simulations by iMachida et al.l ( 120081 ) focus on the outflows/jets at the earliest evolutionary 
phase of star formation, between the prestellar phase and Class phase such as outflows 
from the first adiabatic core and right after 2nd core formation. Theoretical works are able 
to resolve the jet launching points numerically. However, observational studies with current 
telescopes are not able to spatially resolve the region. The rarity of suitable targets (short 
life time at the earliest evolutionary phase) is a persistent problem. 



Takahashi et al 





Menten et al. 


2007 


), is located in the 


3MC- 


Johnstone & Ballv 


1999; 


Matthews et al. 


2005; 



20091 . l201ll ; Hereafter, Takahashi et al. 2011a is denoted as Paper II). 



MMS 6-main is the brightest and the m ost compact submillimeter continu um source in the 



MMS 6 region ( ITakahashi et al.l 120091 ) as well as in the OMC-3 region ( ITakahashi et al. 
201ll ). The sub-arcsecond 850 /zm continuum image (Paper II), spatially resolved a massive 
envelope (0.29 M ) and hot gas (>50 K) in the central 120 AU. Comparisons with core 
models clearly showed the presence of a self-luminous source, which implies that MMS 6- 
main is an intermediate-mass protostellar core at the earliest evolutionary stage. However, 
neither molecular outflow, nor radio jet, nor infrared source has ever been detec ted at the 
source center in the previous studies (IMatthews et al.ll2005t ITakahashi et al.ll2008l . Hereafter, 
Takahashi et al. 2008 is denoted as Paper I). 

In order to search for molecular outflows and study their characteristics, we have per- 
formed sub-arcsecond angular resolution observations toward MMS 6-main in the CO (3-2) 
and HCN (4-3) lines with the SMA. Our observations have achieved a factor of seven im- 
provement in terms of beam size (~50 times better in terms of the beam surface area) than 



the previous molecular line studies by ITakahashi et al. 
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OBSERVATIONS AND DATA REDUCTION 



The observations of the CO (3-2; 345.796 GHz) and HCN (4-3; 354.505 GHz) emission 
were carried out with the SMA (IHo et al.ll2004j) in the extende d configuration on September 
2, 2010. The phase center was set on MMS 6 flchini et aI.lll997L ; R.A. (J2000)=5 /l 35 m 23/48, 
decl. (J2000)=-05°01'32".20). The array configuration provided projected baselines ranging 
from 44-244 kA for CO (3-2) and 45-249 kA for HCN (4-3). Our observations were insensitive 
to str uctures more extended than 3".8 for CO (3-2) and 3".7 for HCN (4-3) at the 10% 
level (IWilner fc Welch! 1 19941 ) . The typical system noise temperatures in DSB mode were 
between 200-350 K at the observed elevations. The receivers have two sidebands with 4 
GHz bandwidth and the CO (3-2) and HCN (4-3) lines were simultaneously observed in the 
lower- and upper-sideband, respectively. We used the configuration that gave 128 channels 
for CO (3-2) and 64 channels for HCN (4-3). The corresponding velocity resolutions were 
0.7 km s _1 for CO (3-2) and 1.4 km s _1 for HCN (4-3). The phase and amplitude calibrator, 
0423-013 (1.5 Jy ), was observed every 18 minutes. Observations of Callisto provided the 
absolute scale for the flux density calibrations. The overall flux uncertainty was estimated to 
be ~20%. The passband across the bandwidth was determined from observations of 3C 454.3 
with a ~30 minute integration. 



The raw data were calibrated using MIR ( iScoville et al.lll993l ). After the calibration and 
the flagging of bad data, final CLEANed images were made using the AIPS task "imager" 
with a natural weighting. The resulting synthesized beam sizes were 0". 82x0". 64 with a 
position angle of -74° for CO (3-2) and -79° for HCN (4-3). The achieved rms noise levels 
were 0.14 Jy beam -1 km s _1 for CO (3-2) and 0.10 Jy beam -1 km s _1 for HCN (4-3), for 
the 0.7 km s -1 bin and 1.4 km s _1 bin. 



3. RESULTS 

Figure 1 shows the zeroth- moment (total intensity), first-moment (velocity field de- 
scribed in the radial outflow velocity), and second- moment maps (velocity dispersion) pro- 
duced in CO (3-2) and HCN (4-3), superposed on the 850 /im continuum emission (Paper 
II). With the SMA extended configuration observations, a very compact molecular outflow 
lobe (~2") associated with MMS 6-main has been successfully imaged and spatially resolved 
for the first time. The outflow has a collimated bipolar structure and the high-velocity 
emission is elongated along the north-south direction with a position angle of -6°. The first 
moment maps show that the redshifted and blueshifted components are clearly separated in 
the north and south directions centered at the MMS 6-main continuum peak. The second 
moment maps show that the line width dramatically increases downstream at the ends of 
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the molecular outflow particularly in the HCN (4-3) emission. This is further demonstrated 
by the position- velocity map presented in Figure 2 in HCN (4-3), where the radial outflow 
velocity reaches its maximum only at ~2" from the exciting source. 

Note that the gas distribution and velocity structure as observed in CO (3-2) and 
HCN (4-3) are slightly different. In Figure 1 (first moment maps), the peak positions of 
the redshifted components are roughly the same in both CO (3-2) and HCN (4-3), while 
the peak positions of the blueshifted component are shifted approximately 1" with respect to 
each other. Moreover, as clearly seen in Figure 2, the CO (3-2) emission traces only a part 
of gas observed in HCN (4-3) with the radial outflow velocity up to 15 km s -1 . No higher 
velocity emission is traced by the CO (3-2) emission. Note that the gas distribution does not 
dramatically change even if we smooth the CO (3-2) data to the same velocity resolution 
as the HCN (4-3) data. It is difficult to explain that the CO (3-2) high-velocity outflow 
is fainter than in HCN (4-3), if the emissions are optically thin, since the CO molecule is 
much more abundant. It is likely that the CO (3-2) is optically thick, and may be sensi- 
tive to the less dense region, perhaps the outflow cavity. If the outflow walls are stratified 
in temperature but not in velocity, self-absorption by the cooler outer layers may suppress 
CO (3-2) relative to the HCN (4-3). In this context, we note that the CO emission region 
is more extensive than in HCN. Multi-transition line observations are necessary to constrain 
detailed physical conditions of the high- velocity component. 

Paper I and Takahashi et al. (2009) did not detect molecular outflows associated with 
MMS 6-main in the CO (1-0) and CO (3-2) observations. This is likely due to the lim- 
ited sensitivity and the beam dilution. The 3 o mass sensitivity limit was 10~ 4 M© (with 
the ASTE 26" beam size; Paper I) and 6.5 xlO" 4 M (with the NMA ~5" beam size; 



Takahashi et al.ll2009i ). Here, with the sub-arcsec resolution (0". 82x0". 64), we detect the 
very compact outflow («2") with the mass sensitivity limit of 3.9 xlO -7 M (3a mass sensi- 
tivity derived from HCN). 



4. Outflow Parameters 

Outflow parameters were calculated from the detected CO (3-2) and HCN (4-3) emis- 
sions. Note that, as we described in Section 2, the maximum detectable size with our SMA 
observations is larger than the outflow lobe size (~2"), so that our observations are not 
significantly affected by missing flux. 
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The inclination angle of the outflow^!, i, is constrained from t he positio n -veloc ity di- 
agram. We use a simple analytical outflow model as described in iLee et al.l (120001 ) . The 
comparison with these models shows that the inclination angle of 45° as denoted by the 
black curves in Figure 2b, can delineate the MMS 6 outflow velocity structure. Therefore, 
the outflow inclination angle of 45° is adopted in the following discussion. 

In order to compare the MMS 6-main outflow properties with those derived from the 
other OMC-2/3 outflows reported in Paper I, the same derivations are adopted. Under the 
assumption of the local thermodynamical equilibrium (LTE) condition and optically thin 
molecular line emissions, we estimate the outflow mass. The excitation temperature of the 
molecular line is adopted as 30K, which is the typical CO (3-2) brightness temperature mea- 
sured in the OMC-2/3 region (Paper I). Note that adopting the excitation temperature of 100 
K modifies the outflow mass estimation slightly (a factor of < 1.5). The abundance of the CO 
and HCN are adopted as XfCOl =lQ- 4 fererking et al.l Il982h and X[HCNl=5xlO~ 7 (values 
derived in the L 1157 outflow by lBachiller fc Perez Gutierrez! 119970. The HCN (4-3) abun- 
danc e can be enhanced by f actors of 10-100 in the shock region (IBachiller fc Perez Gutierrez 



19971 ; iJorgensen et al.ll2004l). In this paper, the maxim um HCN (4-3) abundance derived for 



the L 1157 outflow (IBachiller fc Perez Gutierrez! 119971 ). is adopted. 

The maximum radial outflow velocitj|§ and the outflow lobe size are derived as t>fl ow = 
f max ( bs) x [1/sim] and -Rfl ow = -RobsX [1/cosi], respectively. The dynamical time of the outflow 
is estimated to be td = -Rflow/^flow Y r with the maximum radial outflow velocity. The outflow 
lobe size, and outflow dynamical time, outflow momentum (P = Mco^flow), energy (E = 
M covl ow /2), outflow force (F ob s = P/U), mechanical luminosity (L m = M C o><vl ow /2R), 
and mass loss rate (M out = McoAd) are estimated. The estimated parameters using HCN 
(4-3) and CO (3-2) are listed in Table 1, with/without inclination angle corrections. 



5. Nature of the Extremely Compact Outflow 

In the MMS 6-main outflow case, HCN (4-3) is a better tracer of the morphology and 
velocity structure of the collimated molecular outflow as compared with CO (3-2). This 
may be because the central part of the MMS 6-main envelope (<1500 AU) is extremely 
dense (n~10 9 cm -3 ; Paper II) and hence, the ejected gas is better traced by the higher 
critical density of HCN (4-3), and also because the CO may suffer from self-absorption 



2 Here, the plane of the sky is defined as i=0°. 

3 We correct all velocities by the systemic velocity (see Table 1). 
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effects. Hereafter, we use the outflow properties as estimated from the HCN (4-3) emission. 
Nevertheless, estimated outflow parameters from HCN (4-3) and CO (3-2) do not show 
significant differences. 

The high-velocity gas detected in HCN emission clearly shows the bipolar collimated 
structure centered at MMS 6-main with a projected lobe size of 1700AU (0.0085 pc) for 
the blueshifted component and 1200 AU (0.006 pc) for the redshifted component with an 
inclination angle of 45°. The dynamical timescale of this outflow is estimated to be 33-44 yr 
with the maximum radial outflow velocity of 35-39 km s _1 . (The dynamical time becomes 
a factor of 1.5-2.3 shorter when the intensity weighted mean outflow sizes and mean gas 
velocities are adopted.) The derived outflow mass and outflow size are approximately 2-4 
orders and 1-3 orders of magnitude smaller t han those derived for other molecular outflows 
previously detected in the OMC-2/3 region (lAso et al.ll2000t IWilliams et al.ll2003l , Paper 



I). The molecular outflow detected in MMS 6-main is the smallest and least massi ye bipolar 



outflow that has ever been observed for the intermediate-mass pr otostars (c.f., lAso et al 



200CH : IWilliams et al.ll2003l : IZapata et al.ll2006l : Beltran et all 12008 
outflow force, (9.1-13) xl0~ 5 M km s" 1 yr" 1 



Paper I). In contrast 
and the mass loss rate, 3.5 xl0~ 6 M Q yr -1 



have similar values as compared with the median values/mean value of those derived for 
the other OMC-2/3 outflows (F C o=5.4xlO~ 5 M km s" 1 yr" 1 / 2.7xl0~ 4 M km s~ x yr 
Values are calculated from Paper I). These strongly imply that the outflow associated with 
MMS 6-main is much younger, but otherwise has similar outflow force as compared with the 
other molecular outflows in OMC-2/3. 

The projected outflow lobe size (~1000 AU) is even smaller than the outflows found 
in the protostellar cores in Taurus (6000 AU/18000 AU for the median/mean values; from 
Hogerheijde 1998). The smallest outflow ever reported in the previous studies is an outflow 
associated with a very low-luminosity object (VeLLOs), L 1014L, which has a lobe size of 
~500 AU and outflow mass of <10 -4 M . This source is a proto-brown dwarf candidate 



( IBourke et al.ll2005l ). Detected molecular outflow in MMS 6-main is a factor of three larger 
than this outflow. However, the estimated outflow force in MMS 6-main is more than 
two orders of magnitude larger than that estimated in L 1014L. Furthermore, similarly 
compac t molecular outflows are re ported in young brown dwarf stars; ISO-Oph 102 and 
MH05 flPhan-Bao et al.ll2008[ [jpllh . They have outflow lobe sizes of ~1000 AU and outflow 
masses of ~10 -4 M , which are similar to the outflow in MMS 6-main. However, their 
estimated outflow forces are on the order of 10~ 8 M km s _1 yr -1 , which is four orders 
of magnitude smaller than those estimated in MMS 6-main. As summarized in Figure 3, 
the molecular outflow detected in MMS 6-main is similarly compact and less massive as 
compared to those associated with young brown dwarfs, while the outflow force is much 
higher in MMS 6-main. If we consider that the outflow force is correlated with the mass 
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accretion rate as suggested in lBontemps et al.l (119961 ). this result indicates that MMS 6- main 
has a much higher accretion rate than for (proto) brown dwarf cases. 

The observed outflows are usually charact erized by the jet-driven bow-shock models and 
wind-driven shell models (e.g.. lLee et al.ll2000l ). The jet-driven bow shock model shows that 
molecular outflows consist of the ambient gas interacting with the bow-shock of the jet head, 
while the wind-driven model shows that molecular outflows consist of swept-up gas entrained 
by the wide-angle outflow. The velocity structures observed in HCN (4-3) emission show 
the immediate velocity increment at the jet head up to 25 km s -1 . This clearly suggests the 
signature of the jet-driven bow-shock type outflow. Simil ar velocity structures are reported in 
low-m ass protostellar outflows in the Class phase (e.g. 
2010h . 



Lee et al.ll200d 120071 : iHirano et al. 



The pri mary jets from young stars are likely la unched from accretion disks around pro- 
tostars (e.g.. iPudritz et al.ll2007l ; IShang et al.ll2006l ). Magnetocentrifugal forces can produce 
winds with large opening-angles, while the magnetic pressure gradients (hoop stress induced 
by the generated toroidal magnetic field) event ually dominate bey ond the region around the 
Alfven radius, which provides the collimation ( iPudritz et al.ll2007l ). X-wind models can also 
explain both wide and co llimated winds as an effect due to the density contrast within the 
winds (IShang et al.l 120061 ) . 



Furthermore, recent 3D- MHD simulations pre dict two distinct outflows associated with 
the early evolutionary stages (IMachida et al.l 120081 ). First case is the outflow from the first 
adiabatic core. These outflows have a strong magnetic field and are driven mainly by the 
magneto centrifugal mechanism and are guided by the hour glass-like field lines. These 
outflows show a wide opening-angle low- velocity outflow (~5 km s _1 ). Second case is the 
outflow from the protostar. These outflows have a weak magnetic field. They are driven 
by the magnetic pressure gradient force and guided by straight field lines. They have a 
collimated high- velocity outflow (>30 km s _1 ). 

Our observations do not have enough angular resolution to resolve the outflow launching 
point. Hence, this experiment cannot distinguish between the different launching mecha- 
nisms as proposed by theoretical models. Nevertheless, the nature of the observed outflow in 
MMS 6-main shares similarities with the outflows from protostellar cores, such as collimation, 
bipolarity, high-velocity gas, and bow-shocks. Negative detection of the wide opening-angle 
low-velocity outflow as predicted in the first adiabatic core, is possibly due to the observa- 
tional limits on the mass sensitivity (Section 3) and the detectable size (Section 2). 

The extremely compact and less massive nature of MMS 6-main suggests that it may 
be at the earliest evolutionary phase of the protostar, shortly after its formation. Paper II 
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suggests that the observed density of the submillimeter source, MMS 6-main, is explained by 
the combination of a central heating source and the surrounding warm envelope, supporting 
the presence of a protostellar core. 

6. Conclusion 

An extremely compact molecular outflow was discovered toward an intermediate-mass 
protostellar core, MMS 6-main in the OMC-3 region. The velocity broadening at outflow 
tips, (At>~25 km s _1 ), indicates the bow-shock type of velocity structure. The estimated 
outflow lobe size (~1500 AU), mass(~10 -4 M ), and dynamical time (<100 yr), clearly 
suggest that the outflow is compact, less massive, and young, while the estimated outflow 
force (^10~ 4 M Q km s _1 yr -1 ) is similar to those derived in other OMC-2/3 outflows. These 
results most likely suggest that the outflow is originated from a protostar (2nd core), shortly 
after its formation. These results are consistent with our recent SMA continuum observations 
at 850 jum, which imply that MMS 6-main is an intermediate mass protostellar core. 

We acknowledge the staff at the Submillimeter Array for assistance with observations. 
We thank Naomi Hirano, Masanori Nakamura, and Keiichi Asada for fruitful discussions, 
and our anonymous referee whose suggestions improved this manuscript. 
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Fig. 1. — Color images of the zeroth, first and second moments maps of the HCN (4~3) and 
CO (3-2) emissions. The velocity labels in the first moment maps are described in the radial 
outflow velocity (t> ra d;obs = ^lsr — v sys ). The moment maps were made using the AIPS task 
'momnt' with 5a cut off level (i.e. 0.5 Jy beamT 1 km s _1 for HCN and 7.0 Jy beamT 1 km 
s -\ j or QO) The black contours show the 850 [im continuum image (Paper II). The contour 
levels are 10a, 20a, 30a, 40a, 60a, 100a, 160a, 200a, and 240a with intervals of 20 a 
(la =2.1 mJy beam -1 ). Filled ellipses at the bottom-left corners show the synthesized beams 
of the line data (Section 2) and continuum data (0". 35x0". 30 P. A. =32°). 
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Fig. 2. — (a) HCN(4~3) and CO (3-2) position-velocity maps denoted by black and purple 
contours, respectively. The velocity labels are described in the radial observed velocity (^lsr/ 
These position-velocity maps were made along the major axis (P. A. =-6° ) of the molecular 
outflow. Dashed lines along the horizontal and vertical directions show the position of the 
850 fim continuum peak (Paper II) and the systemic velocity of MMS 6 as derived from 
the H 13 CO + (1-0) emission (Takahashi et al. 2009), respectively. The resolution of each 
molecular line in the position-velocity diagrams are denoted in the bottom right corner, (b) 
HCN (4~3) position-velocity map as shown in figure (a) overlaid with the wide-angle wind 
model curve produced by i=45°, C=0.88 arcsec~ l , vq=7.0 km s _1 (refer Lee et al. 2001 for 
definitions of C and v ). This model is used in order to constrain the inclination angle of 
the molecular outflow. 
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Fig. 3. — a: Outflow mass plotted as a function of outflow size, b: Outflow force plotted as 
a function of outflow size, c: Outflow force plotted as a function of bolometric luminosity. 
Connected lines of data points show the uncertainties of bolometric luminosities measured 
by Paper I. Black and blue arrows show the lower and upper limits of the bolometric lu- 
minosities measure by Paper II and Chini et al. 1997, respectively. Blue filled diamond 
denotes outflow from MMS 6-main. Open circles, filled circles, open squares, crosses, open 
triangle s, and hatched squar es show outflow results reported towards intermediate-mass pro- 
tostars iBeltrdn et al\ \200a, and their references ), OMC-2/3 outflows (in the blueshifted 
and redshifted compone nts; Takahashi et al. 200s]) , Taurus outflows (in the blueshifted and 



redshifted components; \Hogerheijde et al 



1 998 ) , and outflows associated with proto - brown 



dwarf and first ad iabatic core candidates iBourke et all \200a ; \Phan-Bao et all \2008\ . \20lA ; 



Pineda et al\\20lA) . respectively. Data points from Taurus and OMC-2/3 surveys are sep- 
arately plotted in redshifted and blueshifted components. Adding both components should 
increase the values by a factor of two as compared to the current plot. 
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Table I. Outflow Parameters Estimated from HCN (4-3) and CO (3-2) 



.Parameters 




.Diue 


Lobe 


irteo 


Lobe 






£ — *±o 


TT 

U ncorrcctcci 




TT 

U ncorrcctcci 






HCN (4 


-3) 






Mass (Mq) 




8.3e-05 


8.3e-05 


5.2e-05 


5.2e-05 


Maximum velocity (km s _1 


) a 


39 (20) b 


27 (14) b 


35 (18) b 


25 (13) b 


Size (AU) 




1700 (1100) c 


1200 (780) c 


1200 (740) c 


830 (520) c 


Dynamical Time (yr) 




44 (56) d 


44 (56) d 


33 (40) c 


33 (40) d 


Momentum (Mq km s _1 ) 




3.2e-04 


2.3c-04 


1.8c-04 


1.3e-04 


Kinetic energy (Mq km 2 s" 




6.2e-03 


3.1c-03 


3.2c-03 


1.6e-03 


Outflow force (Mq km s~ 1 


yr- 1 ) 


7.3e-05 


5.2e-05 


5.5e-05 


3.9c-05 


Mechanical luminosity (Lq 


) 


4.9e-03 


2.4c-03 


3.4e-03 


1.7e-03 


Mass Loss rate (Mq yr -1 ) 




1.9e-06 


1.9e-06 


1.6e-06 


1.6c-06 






CO (3- 


2) 






Mass (M Q ) 




1.9e-04 


1.9c-04 


1.4e-04 


1.4e-04 


Maximum velocity (km s _1 


) a 


21 (ll) b 


15 (7.5) b 


11 (5.5) b 


7.7 (3.9) b 


Size (AU) 




1300 (570) c 


950 (400) c 


1100 (640) c 


770 (450) c 


Dynamical Time (yr) 




64 (53) d 


64 (53) d 


100 (115) d 


100 (115) d 


Momentum (Mq km s -1 ) 




4.1e-04 


2.9c-04 


1.5c-04 


l.lc-04 


Kinetic energy (Mq km 2 s~ 


- 2 ) 


4.3e-03 


2.1c-03 


8.2c-04 


4.1c-04 


Outflow force (Mq km s~ 1 


yr- 1 ) 


6.3e-05 


4.5e-05 


1.5e-06 


l.lc-05 


Mechanical luminosity (Lq 


) 


2.3e-03 


1.2e-03 


2.9e-04 


1.4e-04 


Mass Loss rate (Mq yr -1 ) 




3.0c-06 


3.0e-06 


1.4c-06 


1.4e-06 



a The maximum radial outflow velocities (f Q bs — I^LSR. — ^sysl) measured from the channel maps. 
Here, v sys = 10.3 km s- 1 is adopted (Takahshi et al. 2009) 

b Mcan radial outflow velocities. 

c Distance to the peak intensity positions from the driving source measured from the zeroth 
moment maps 

d Dynamical times derived with the mean gas velocities and the outflow intensity peak positions. 



